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ff) (57) Abstract: A filter settings generation operation includes sampling a communication channel to produce a sampled signal. The 
V) sampled signal is spectrally characterized across a frequency band of interest to produce a spectral characterization of the sampled 
0\ signal. This spectra] characterization may not include a signal of interest. The spectral characterization is then modified to produc 
_ a modified spectral characterization. Filter settings are then generated based upon the modified spectral characterization. Finalh 



a- settings are then generated based upon the modified spectral characterization. Finally, 
^ M e communication channel is filtered using the filter settings when the signal of interest is present on the communication channel. 

In modifying the spectral characterization, pluralities of spectral characteristics of the spectra] characterization are independently 
£jj modified to produce the modified spectral characterization. Modifications to the spectral characterization may be performed in 
the frequency domain and/or the time domain. One particular spectra] modification that is performed is raising of the noise floor 
Q of the spectral characterization to meet a budgeted signal-to-noise ratio. Other spectral modifications include modifying spectral 
^ components corresponding to an extended interfering signal. In modifying these spectral characterizations, spectral components 
^ corresponding to a plurality of expected interfering signals may be modified. 
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TITLE: SYSTEM AND METHOD FOR CANCELING INTERFERENCE IN A 
COMMUNICATION SYSTEM 

SPECIFICATION 

5 1. FIELD OF THE INVENTION 

The present invention relates generally to communication systems, and more 
particularly to interference mitigation in communication systems. 
2. BACKGROUND OF THE INVENTION 

The structure and operation of high-speed data communication systems are generally 
1 0 known. Such high-speed data communication systems employ various media and/or wireless 
links to support the transmission of high-speed data communications. Particular 
embodiments of high-speed communication systems include, for example, cable modem 
systems, home networking systems, wired local area networks, wired wide area networks, 
wireless local area networks, satellite networks, etc. Each of these high-speed data 

1 5 communication systems has some unique operational characteristics. Further, some of these 
high-speed data communication systems share similar operational drawbacks. Home 
networking systems and cable modem systems, for example, are both subject to interfering 
signals that are coupled on media that carry communication signals. 

Cable modem systems, and more generally, cable telecommunication systems include 

20 set-top boxes and residential gateways that are in combination capable of currently providing 
data rates as high as 56 Mbps, and are thus suitable for high-speed file transfer, video 
teleconferencing and pay-per-view television. These cable telecommunication systems may 
simultaneously provide high-speed Internet access, digital television (such as pay-per-view) 
and digital telephony. Cable modems are used in a shared access environment in which 

25 subscribers compete for bandwidth that is supported by shared coaxial cables. During normal 
operations, sufficient bandwidth is available across the shared coaxial cables to service a 
large number of subscribers, with each subscriber being serviced at a high data rate. Thus, 
during normal operations, each subscriber is provided a data rate that is sufficient to service 
uninterrupted video teleconferencing, pay-per-view television, and other high bandwidth 

30 services. 

Intermittent, narrowband interfering signals may, from time to time, interfere with 
wideband communication signals, e.g., upstream Data-Over-Cable Service Interface 
Specification (DOCSIS) transmissions ("desired signals"). These intermittent narrowband 
interfering signals unintentionally couple to the shared coaxial cables via deficiencies in 
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shielding and/or other coupling paths. With these interfering signals present, the data rate 
that is supportable on the coaxial cables is reduced, hi some cases, depending upon the 
strength and band of the interfering signals, the supportable bandwidth is reduced by a 
significant level. 

5 Conventionally, when an interfering signal is present, an adaptive cancellation filter is 

employed by each cable modem receiver to cancel the interfering signal by adoptively 
placing a filtering notch or null at the frequency of the interfering signal. When the 
interfering signal becomes absent, the conventional adaptive cancellation filter continues to 
adapt and removes the filtering notch. If the mterfering signal reappears, the adaptive 

1 0 cancellation filter again adapts to null the interfering signal. Thus, when the interfering 

signal first reappears, the cancellation filter cannot fully compensate for the interfering signal. 
Because many interfering signals are intermittent, the presence of these intermittent signals 
reduces the bandwidth that is supportable upon the coaxial media during the time period 
required for the cancellation filter to adapt. Further, because the interfering signals 

1 5 oftentimes vary in strength while present, the adaptive cancellation filter most often times 
does not fully remove the interfering signal. 

Overlapping adjacent channel signals cause another source of interference for the 
desired signal because they often produce interfering signals in the band of the desired signal. 
For example, a TDMA signal that resides in an adjacent channel and that turns on and off 

20 may have side lobes that overlap and interfere with the desired signal. When the interfering 
signal is present, the. conventional adaptive cancellation filter places a notch or null at the 
frequency band of the interfering signal. When the interfering signal is absent, the 
cancellation filter adapts to removes the notch. The precise amount of interference in the 
desired signal caused by the adjacent channel signals may vary with data content in the 

25 adjacent channel. 

Thus, in both the case of the narrowband interferer and the adjacent channel 
interferer, the interfering signal(s) varies over time. For this reason, an optimal or near- 
optimal solution may be found only for the average interfering strength of the interfering 
signal(s), but not for the peak(s) of the interfering signal(s). In many operational conditions, 

30 typical fluctuations in the strength of interfering signals cause conventional cancellation 
filters to provide insufficient cancellation. Resultantly, overall bandwidth that could be 
provided by the supporting communication system on the particular shared media is 
significantly reduced. 
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Therefore, there is a need in the art for a filtering system and associated operations 
that cancel interfering signals so that throughput is maximized. 

SUMMARY OF THE INVENTION 

5 The present invention specifically addresses and alleviates the above-mentioned 

deficiencies associated with the prior art to efficiently filter a communication channel to 
remove narrow band interfering signals. According to the present invention, a 
communication channel is sampled to produce a sampled signal. Then, the sampled signal is 
spectrally characterized across a frequency band of interest to produce a spectral 

10 characterization of the sampled signal. This spectral characterization may not include the 
signal of interest. The spectral characterization is then modified to produce a modified 
spectral characterization. Filter settings are then generated based upon the modified spectral 
characterization. Finally, the communication channel is filtered using the filter settings when 
the signal of interest is present on the communication channel. 

15 Various operations are employed to modify the spectral characterization to produce 

the modified spectral characterization. In most cases, a plurality of spectral characteristics of 
the spectral characterization is independently modified to produce the modified spectral 
characterization. Modifications to the spectral characterization may be performed in the 
frequency domain and/or the time domain. One particular spectral modification that is 

20 performed is raising of the noise floor of the spectral characterization to meet a budgeted 
signal-to-noise ratio. 

Other spectral modifications include modifying spectral components corresponding to 
an expected interfering signal. In modifying these spectral characterizations, spectral 
components corresponding to a plurality of expected interfering signals may be modified. In 

25 obtaining information to perform this spectral modification, spectral components of prior 
sampled signals may be employed. For example, an interfering signal may be intermittent so 
that it is not present in the current sampled signal but was present in prior sampled signals. 
Since overall system performance may be enhanced when the filter compensates for this 
interfering signal even when it is not present, in modifying the spectral characterization of the 

30 sampled signal to produce the modified spectral characterization, the presence of an 
interfering signal in a prior sampled signal may be weighted more heavily than the absence of 
the mterfering signal in the current sampled signal. 



3 
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The communication channel may be sampled with or without the presence of the 
signal of interest. However, when the signal of interest is present on the communication 
channel during the sampling interval, the signal of interest must be removed from the 
sampled signal. Further, when the total spectral density of the sampled signal exceeds a 
5 threshold value, the sampled signal may be discarded as invalid. 

Because the filtering operations of the present invention adapt to time varying 
interfering signals in an adaptive fashion, more efficient filtering operations are performed. 
For example, by raising the noise floor of the spectral characterization, budgeted SNR design 
constraints may be met with minimal filter tap magnitudes. Further, by enabling the 
10 introduction of known interferers into the modified spectral characterization, the filtering 
operations will compensate for mterfering signals immediately when they appear, not after a 
learning period of the prior filtering operations. Because of these benefits, among others, 
greater system throughput is realized. 

Other features and advantages of the present invention will become apparent from the 
15 following detailed description of the invention made with reference to the accompanying 
drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

These and other features, aspects and advantages of the present invention will be more 
20 fully understood when considered with respect to the following detailed description, 
appended claims and accompanying drawings wherein: 

FIG. 1 is a simplified block diagram illustrating a first embodiment of a narrow band 
interference cancellation filter (ICF) used in accordance with the principles of the present 
invention; 

25 FIG. 2 is a simplified block diagram illustrating a set of operations employed to 

generate the filter settings of FIG. 1; 

FIG. 3 is a block diagram illustrating operations employed to generate filter settings 
corresponding to a first embodiment of FIG. 2; 

FIG. 4 is a block diagram illustrating operations employed to generate filter settings 
30 corresponding to a second embodiment of FIG. 2; 

FIG. 5 is a graph illustrating the spectral characteristics of a sampled signal according 
to the present invention; 



4 



WO 02/101953 



PCT/US02/18521 



FIG. 6 is a graph illustrating the modified spectral characteristics of a sampled signal 
according to the present invention; 

FIG. 7A is a simplified block diagram illustrating an embodiment of a narrow band 
interference cancellation filter used in accordance with the principles of the present invention; 
5 FIG. 7B is a block diagram of a filter weight computation algorithm for the ICF of 

FIG.7A; 

FIG. 8 is a logic diagram illustrating operations performed corresponding to the filter 
weight computation algorithm of FIG. 7B; 

FIG. 9 is a table of modes used in connection with a spectral mask incorporated in the 
10 algorithm of FIG. 7B; 

FIG. 10 is a diagram illustrating a non-linear function used in connection with the 
tracker of the filter weight computation algorithm of FIG. 7B; 

FIG. 11 is a flow diagram depicting a noise floor setting algorithm for the filter 
weight computation algorithm of FIG. 7B; 
15 FIG. 12 is a flow diagram depicting a modified noise floor setting algorithm for the 

filter weight computation algorithm of FIG. 7B; and 

FIG. 13 is a system diagram illustrating the manner in which the present invention 
may be employed in a cable telecommunication system servicing hi$i speed data 
communications. 

20 

DETAILED DESCRIPTION OF THE INVENTION 

FIG. 1 is a simplified block diagram illustrating a first embodiment of a narrow band 
interference cancellation filter (ICF) used in accordance with the principles of the present 
invention. According to the present invention, incoming unfiltered commumcations 102 are 

25 received on a communication channel and filtered using filter 104 to produce filtered 
communications 110. Filter settings 108, generated by filter settings generation block 106, 
when applied to filter 104 remove narrow band interfering signals from the communication 
channel to produce the filtered communications 110. Thus, during normal operations, the 
filtered communications 1 10 are provided to the receiver section of a device in which the ICF 

30 resides. Because of the filtering operations of the present invention, the receiver is better able 
to operate on a signal of interest (desired signal) contained in the filtered communications 
1 10 to remove data contained therein. 
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FIG. 2 is a simplified block diagram illustrating a set of operations employed to 
generate the filter settings 108 of FIG. 1. In operation, the filter settings generation block 106 
first samples the communication channel 200 to produce a sampled signal 204. The signal 
sampler 202 of FIG. 2 performs sampling operations. In the particular embodiment of 
5 operation illustrated, the communication channel 200 is sampled when a signal of interest 
(desired signal) is not present on the communication channel 200. However, in other 
operational embodiments, the signal of interest is present on the communication channel 200 
during sampling but all or a portion of the signal of interest is removed from the sampled 
signal 204 by the signal sampler 202, taking advantage of a known preamble, for example. 

10 The sampled signal 204 is then spectrally characterized across a frequency band of 

interest by spectral characterization block 206 to produce a spectral characterization 208 of 
the sampled signal 204. Because the sampled signal 204 did not include the signal of interest, 
the spectral characterization 208 also does not include the signal of interest 212. Thus, the 
sampled signal 204 includes noise that is present during the sampling period and includes 

15 interfering signal(s) if present. The filter settings generation block 106 operates periodically 
to generate new filter settings 108. Thus, during some sample periods, interfering signals are 
present while during other sample periods, interfering signals are not present. Therefore, for 
any given operation of the filter settings generation block, the spectral characteristics 208 
may or may not include interfering signals. However, the spectral characteristics 208 will 

20 always include noise at some level. 

Next, the spectral characteristics 208 of the sampled signal 204 are modified to 
produce a modified spectral characterization 212 by modify spectral characterization block 
210. The modifications performed by modify spectral characterization block 210 to produce 
the modified spectral characterization 212 are numerous and varied. These operations 

25 include raising the noise floor of the spectral characterization 208, introducing spectral 
characteristics into the spectral characterization 208 that relate to previously present 
interfering signals, introducing spectral characteristics of expected interfering signals into the 
spectral characterization 208 that are expected but that have not recently been present (or at 
all been present), and other preemptive modifications. 

30 A filter settings generation block 214 men generates filter settings 1 08 based upon the 

modified spectral characterization 212. Finally, the filter settings 108 are applied to the filter 
104 for subsequent filtering of the communication channel when the signal of interest is 
present on the communication channel. 
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FIG. 3 is a block diagram illustrating operations employed to generate filter settings 
corresponding to a first embodiment of FIG. 2. As shown in FIG. 2, the communication 
channel is sampled using signal sampler 302 to produce a sampled signal 304 in the time 
domain. The signal sampler 302 produces the sampled signal 304 so that a signal of interest 
5 is not present (or substantially attenuated). The sampled signal 304 is a digital representation 
of the analog signal that is present on the communication channel 300 during the sampling 
interval. 

A Fast Fourier Transform (FFT) is then performed upon the sampled signal 304 by 
spectral characterization block 306. A Harming or other window is used to remove 

10 sensitivity to frequency location. The FFT operation (as well as most/all other of the 
operations upon the digital data) is performed using a Digital Signal Processor (DSP) or 
another computing device resident in the device performing the operations of the present 
invention. The spectral characterization block 306 perfonriing the FFT operation therefore 
produces a spectral characterization 308 of the sampled signal. Because no signal of interest 

15 (desired signal) was present in the sampled signal 304, no spectral components of the signal 
of interest (desired signal) are present in the spectral characterization 308. 

As is illustrated, the spectral characterization block 306 of FIG. 3 receives input in the 
time domain but produces output in the frequency domain. The concepts of transformation 
from the time domain to the frequency domain and the transformation from the frequency 

20 domain to the time domain are well known and will be described herein only as they relate to 
the present invention. Modification of the spectral components of the sampled signal may be 
performed in the time domain and/or the frequency domain. In the embodiment of FIG. 3, 
modification of spectral components is first performed in the frequency domain by modify 
spectral characterization block 3 10 to produce a modified spectral characterization 3 12. Such 

25 modification of the spectral components by the modify spectral characterization block 310 
may include raising the noise floor of the spectral characterization 308, introducing spectral 
characteristics that relate to previously present interfering signals, and introducing spectral 
characteristics of expected interfering signals, among others. 

The modified spectral characterization 312 is then converted back to the time domain 

30 as modified spectral characterization 316 via Inverse Fast Fourier Transform (IFFT). In the 
time domain, the modified spectral characterization 316 is contained in an autocorrelation 
matrix, which is applied to a filter settings generation block 318 that produces the filter 
settings 320. Note that during the generation of the filter settings, the modified spectral 
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characterization 316 of the sampled signal may be further modified. The modified spectral 
characterization 316 is further modified by modifying the components of the corresponding 
autocorrelation matrix in the time domain. These further modifications may be required to 
generate filter settings that meet a set of filter settings criteria. Such further modification (as 
5 will be described further with reference to FIGs. 7B and 8) may include further raising the 
noise floor of the spectrum corresponding to the modified autocorrelation characterization 
3 1 6 and using prior filter settings in order to generate a satisfactory set of filter settings. 

FIG. 4 is a block diagram illustrating operations employed to generate filter settings 
corresponding to a second embodiment of FIG. 2. As contrasted to the operations of FIG. 3, 

10 the operations of FIG. 4 are performed fully in the time domain. As shown in FIG. 3, the 
communication channel is sampled using signal sampler 412 to produce a sampled signal 414 
in the time domain that either does not contain the signal of interest or from which the signal 
of interest was substantially attenuated. This sampled signal 414 is then presented to 
autocorrelation block 416 that determines the autocorrelation function of the sampled signal 

15 414. The result of this operation is the generation of a specttal characterization 418 in the 
time domain of the sampled signal 414 in the form of an autocorrelation matrix. 

The modify autocorrelation characterization block 420 operates on the spectral 
characterization 418 of the sampled signal 414 in the time domain by altering the coefficients 
of the autocorrelation matrix. The altered autocorrelation matrix therefore represents the 

20 modified spectral characterization 422 that is in the time domain. Modification of the 
spectral components of the spectral characterization 418 by the modify autocorrelation 
characterization block 420 may include raising the noise floor of the spectral characterization 
418, introducing spectral characteristics that relate to previously present interfering signals, 
and introducing spectral characteristics of expected interfering signals that have not yet been 

25 present, among others. The noise floor of the modified spectral characterization 418 of the 
sampled signal 414 may be changed by altering the diagonal components of the 
autocorrelation matrix. Other spectral characteristics of the spectral characterization 418 of 
the sampled signal 414 may be changed by altering the off-diagonal components of the 
autocorrelation matrix. 

30 The modified spectral characterization 422 is applied to a filter settings generation 

block 424 that produces the filter settings 426. Note that during the generation of the filter 
settings, the modified spectral characterization 422 of the sampled signal may be further 
modified. Such further modification may be required to generate filter settings that meet a set 
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of filter settings criteria. Such further modification (as will be described further with 
reference to FIGs. 7B and 8) may include further raising the noise floor of the modified 
spectral characterization 422 and using prior filter settings in order to generate a satisfactory 
set of filter settings. Spectral characteristics of the modified spectral characterization 422 are 
5 further altered in the filter settings generation block 424 by modifying the components of the 
corresponding autocorrelation matrix. 

FIG. 5 is a graph illustrating the spectral characteristics of a sampled signal according 
to the present invention. The representation of the sampled signal of FIG. 5 is in the 
frequency domain, e.g., spectral characterization 308 of FIG. 3. In one embodiment of the 

10 present invention, the spectral characterization will include 256 unique "frequency bins" 
within a frequency band of interest. Each frequency bin corresponds to a unique frequency, 
includes a magnitude component and may include an angular component. In a particular 
implementation of the present invention, the represented magnitude represents the power of 
the sampled signal, e.g., signal magnitude squared, without an angular component. Further, 

15 in the particular embodiment, the represented magnitude is expressed in a logarithmic scale, 
e.g., log base 2, for each frequency bin. 

Thus, in the graph of FIG. 5, the frequency bins include components for 1, 2, . . . , N, 
where N = 256. The spectral characterization of the sampled signal includes a spectral 
component for each frequency bin, denoted as Sfl, . . ., Se*. Because no signal of interest 

20 (desired signal) is present in the spectral characterization of FIG 5, each spectral component 
corresponds to either noise or a combination of noise and a sampled narrowband interferer. 
For example, spectral components Sfl and Sf4 correspond to a particular narrowband 
interferer that was present on the communication channel during the sampling interval. Other 
spectral components corresponding to interfering signals that are present during the sampling 

25 interval may also be present but are not shown in the example of FIG. 5 . 

Also illustrated in FIG. 5 is a representation of a noise floor for the spectral 
characterization. The level of the noise floor may be calculated as an average of spectral 
components that do not exceed a threshold, an average of all spectral components, a weighted 
average of all or a portion of the spectral components, or in another manner that accurately 

30 represents the noise floor for the sampled signal. 

Also illustrated in FIG. 5 is the concept of a budgeted signal to noise ratio (SNR) for 
the communication channel. As is generally known, in the design of a receiver, a minimum 
SNR is desired for satisfactory demodulation operations to be performed. This minimum 
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SNR is referred to herein as the budgeted SNR. Based upon this budgeted SNR and an 
estimate of the signal strength of the signal of interest (desired signal) that will be received by 
the receiver, an estimated level of an acceptable noise floor is determined. As is illustrated in 
FIG. 5, the sampled noise floor lies below the difference between the magnitude of the signal 
5 of interest (desired signal) and the budgeted SNR. Thus, in the example of FIG. 5, the 
receiver could operate within its design criteria even if the noise floor of the communication 
was higher than the sampled noise floor. 

FIG. 6 is a graph illustrating the modified spectral characteristics of a sampled signal 
according to the present invention. As compared to the spectral characterization of FIG. 5, 

10 the modified spectral characterization includes a raised noise floor and an added narrowband 
interferer at frequency bin f7. In raising the noise floor, the magnitude of some, most, or all 
of the spectral components is raised until the raised noise floor corresponds to the difference 
between the expected magnitude of the signal of interest (desired signal) and the budgeted 
SNR. In adding a spectral component of the added narrowband interferer, the magnitude of 

15 one or more spectral components corresponding to the added narrowband interferer is 
increased to a desired magnitude. The modified spectral components are then used to 
generate the filter settings. 

The spectral component modifications illustrated with reference to FIGs. 5 and 6 are 
performed in the frequency domain. However, similar operations may also be performed in 

20 the time domain by operating upon the autocorrelation matrix, as was previously discussed. 
Further, according to one particular embodiment of the present invention (that will be 
discussed with reference to FIGs. 7A, 7B, and 8) the spectral characterization is first 
modified in the frequency domain and then, depending upon whether such modification is 
required to meet filter design criteria, additional spectral modifications are made in the time 

25 domain. 

FIG. 7A is a simplified block diagram illustrating an embodiment of a narrow band 
interference cancellation filter used in accordance with the principles of the present invention. 
In FIG. 7A, a narrow band interference cancellation filter (ICF) 10 is operative to null out 
narrowband interference such as ingress and/or adjacent channel interference from an input 
30 signal. In the embodiment, ICF filter 10 generates a set of adaptive, predictive filter tap 
weights, based on a weight computation algorithm. In the embodiment, the ICF tap weights 
are computed during idle slots when interference is present, but no desired signal is present. 



10 



WO 02/101953 



PCT/US02/18521 



ICF filter 10 includes a delay element 12 to provide a unit delay of one sample, an 
adjustable finite impulse response (FIR) filter 14, and a tap weight computation circuit 16. 
An input signal having narrow band interference such as ingress or adjacent channel 
interference is received by ICF filter 10. Hie input signal is introduced to delay element 12, 
5 tap weight computation circuit 16, and a summing junction 18. Circuit 16 computes tap 
weights, which adjust the pass band characteristics of FIR filter 14 to null out, i.e., cancel, 
narrow band interference appearing in the input signal. For example, if a noise spike or an 
intermittent carrier from another channel appears in the signal, FIR filter 14 would be 
adjusted to have a notch at the frequency of the interference appearing in the input signal. 

10 The output of delay element 12 is coupled to FIR filter 14. FIR filter 14 filters the input 
signal based on the tap weights. In a typical embodiment of the invention, FIR filter 14 
might have 16 taps. The output of FIR filter 14 is combined with the input signal in summing 
junction 18. The output of summing junction 18 is thus equal to the input signal with the 
interference attenuated (or nulled). 

15 The signal distortion caused by the ICF 10 may be compensated by a decision 

feedback equalizer (DFE) 15 in combination with a summing node 19. The filter taps 
generated by the compute coefficients block 16 are also provided to the DFE 15. The DFE 
15 receives the filter taps and, based upon the filter taps, generates a compensation output 
that is summed with the input signal with ingress attenuated at sununing node 19. The output 

20 of the summing node 19 is then received by a decision block 17, which also receives input 
from the DFE 15. The decision block 17, based upon the input from the summing node 19 
and the DFE 15, produces received data. 

In FIG. 7B, the functional elements of the algorithm represented as tap weight 
computation circuit 16 of FIG. 7 A are shown in more detail. The input signal from FIG. 7 A, 

25 which includes narrow band interference such as ingress and adjacent channel disturbances, 
is processed as shown in FIG. 7B. First, the narrow band interference is separated from the 
signal. In the embodiment, a block of 256 complex samples is taken from the communication 
channel during idle Time Division Multiple Access (TDMA) slots, thereby representing the 
narrowband interference. For example, as illustrated in co-pending application Serial No. 

30 09/710,238, filed on November 9, 2000, the MAC at a cable modem termination system 
could create idle slots that insure there are time intervals in which no cable modems are 
fransrnitting an upstream signal. The samples are taken during these idle slots. In another 
embodiment, the interference is separated from the signal by estimating the signal (or a 
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known portion thereof, such as a preamble) and subtracting the estimated signal from the 
input signal, leaving only the narrow band interference. In yet another embodiment, periods 
of no signal on the communication channel are sensed by a threshold detector or SNR 
detector, and the samples are only taken during such periods. In any case, time intervals 
5 containing narrow band interference without signal are produced for sampling and processing 
by circuit 16. 

Preferably, the input samples are converted to the frequency domain using a fast 
Fourier transform (FFT) block 22, which takes a complex fast Fourier transform of the input 
samples. Working in the frequency domain facilitates a number of the operations described 

10 below such as verifying the narrow band content, tracking the changes in the interference 
with different time constants, and frequency masking. However, the tap values could be 
computed in the time domain using autocorrelation properties in another embodiment (as 
described in FIG. 4). In the embodiment of FIG. 7B, a 256-point windowed FFT of the input 
samples is taken to produce a representation of the sampled spectrum without signal, i.e. the 

15 spectrum of the narrow band noise only. Typically, as many as 100 such blocks of samples 
are captured during a time interval T (typically, T=l second), so as to permit rapid tracking of 
time varying narrow band interference. 

The sampled spectrum without signal comprises 256 frequency bins each having 
complex values that represent the energy in the frequency band of the bin. These frequency 

20 bins were illustrated generally with reference to FIGs. 5 and 6. Preferably, the complex 
values of each bin are converted to the absolute value of the result that is transformed to the 
logarithm of the absolute value, thus providing an estimate of the log power spectral density 
(PSD) of the interference. 

To verify the narrow band content of the sampled energy, the output of block 22 is 

25 summed over all the frequency bins comprising the log spectrum, and this sum is compared 
to a threshold value, as represented by a block 30. If the sum exceeds the threshold value, 
then that spectrum is discarded, and no further processing is performed during that particular 
time slot. The spectrum is discarded if the threshold is exceeded because that indicates that a 
large amount of the energy of the spectrum comprises strong wide band energy, or other 

30 content that is inconsistent with narrowband interference. If such a spectrum were included 
in the computations made to generate the tap weights, instead of being discarded, it could 
result in corrupting the narrow band interference spectrum estimate, and therefore the tap 
weights. 
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If the sum passes the threshold test, the spectrum is then introduced to a maximum- 
hold tracking filter represented by a block 32, which tracks the changes in the interference 
with different time constants. This tracking filter maintains an upper bound (with leak) on 
each frequency bin, over multiple input spectra, i.e., multiple blocks of input samples. 
5 Preferably, the tracking filter implements a fast-attack, slow-decay function, or a leaky 
maximum-hold function, depending on the parameters of the error gain function f(x) applied 
by filter 32. As illustrated in FIG. 10, the error gain function f(x) is preferably a nonlinear 
function with a high gain for large positive errors, thereby causing it to rapidly increase when 
a new narrow band interference source appears, and a very low gain for large negative errors, 

10 so that the spectrum estimate remembers an interfering signal for a period of time after the 
narrow band interference disappears. Therefore, in the case of intermittent narrow band 
interference at a given frequency, the filter does not require significant retraining every time 
the narrow band interference reappears. In summary, the tracking filter produces an average 
of the input spectra taken over different periods of time, depending upon the error gain 

15 function. 

Preferably, different gain characteristics similar to those shown in FIG. 10 are .applied 
to the frequency bins individually or in groups. For example, eight different characteristics 
could be used and three control bits could be assigned to each bin in order to select the gain 
characteristic for that bin. Depending upon the particular interference pattern of the 

20 communication channel, the control bits are selected to combat the interference. The output 
of the tracking filter block 32 is then applied to a frequency mask represented by a block 44, 
where a user-defined spectrum mask may be applied to cancel selected frequency bins of the 
filtered spectrum leaving the tracking filter. 

Each of the frequency bins of the mask has associated with it two mode bits, 

25 designating how the mask and the spectrum are to be combined. (See FIG. 9.) The user- 
defined mask may indicate the presence of overlapping, adjacent-channel energy, or a known, 
narrowband interferes, or the like. The spectrum of the permanent energy spectrum of the 
channel can be monitored either by a human observer or automatically by spectrum analysis 
and intelligent computer algorithms, to identify for example adjacent overlapping channel 

30 energy or known interference sources; the mask cancels these unwanted energy components 
in a flexible way based on the selected mode. The user may select from a number of different 
modes of combining the inputs to block 44 on a bin by bin basis. For example, the input 
spectrum could be used en toto and the mask disregarded; the mask could be used en toto and 
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the input spectrum disregarded; the larger of the input spectrum or the mask could be used or 
the sum of the input spectrum and the mask in dB could be used. Each combining method is 
selected independently from bin to bin. 

The noise floor of the spectrum is also adjusted according to the present invention in 
5 mask block 44. Such noise floor alteration was described with reference to FIGs. 5 and 6 
and, while performed by block 44 of FIG. 7B, may be performed within or adjacent to tracker 
block 32. The concept of an SNR goal is used to determine how high the noise floor should 
be in a robust design. For example, for a communications system that uses uncoded 16-QAM 
modulation, a goal of 20 dB is required for good bit error rate (BER) performance. If the 
1 0 thermal noise floor is very low, for example if 30 dB energy per symbol to noise density ratio 
(Es/No) is available, then it is not necessary to null out the narrow band interference all the 
way down to the white noise floor. For example, a 0 dBc interferer only needs to be nulled 
by a little over 20 dB, not by 30 dB. In this example, one would raise the noise floor by 
approximately 10 dB. 

1 5 The noise floor of the spectrum is raised to the SNR goal for the following reasons: 

• It tends to make the numerical computations of the filter coefficients more robust. 

• It helps whiten the noise floor, de-emphasizing small peaks that may have been caused by 
the tracking filter. 

• It tends to reduce the size of the ICF tap weights, which is desirable from the standpoints 
20 of demodulator tracking and DFE error propagation. 

In adjusting the noise floor via mask block 44 operations, the spectrum is processed to 
find those frequency bins that correspond to interference peaks, and those that comprise a 
noise floor. One way to perform this processing is by simple comparison with a threshold. 
Another way (preferred) is by sorting the spectrum into peaks and floor on a bin by bin basis. 

25 The floor bins are averaged in dB and adjusted to produce the desired SNR goal. The 
spectrum intensity in these noise floor, bins (those below the threshold) is then further 
adjusted to the point where the budgeted SNR goal of the communication system is just met. 
It is assumed in this discussion that the spectrum being processed contains higher SNR than 
the SNR goal; if not, insufficient margin exists and the noise floor may not be beneficially 

30 raised. 

After application of the mask and raising of the noise floor, the spectrum is converted 
to its anti-log and the inverse FFT function is performed on the output of block 44 to return to 
the time domain as represented by an autocorrelation block 46. If n is the number of taps in 
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FIR filter 14, the first n time bins of the IFFT at the output of block 46 are used to produce 
the coefficients of the taps of filter 14. After computation of the tap weights and comparison 
with constraints, it is determined if further modifications to the noise floor should be made; 
this is addressed in Figures 8, 11, and 12. If further changes in the noise floor are needed, 
5 then, as represented by R(0) adjust block 48, these changes may be accomplished via 
adjustment of the R(0) term of the autocorrelation function, which is output by block 46. 

The noise floor from block 48 and the first n-1 time bins of the IFFT are supplied to a 
prediction algorithm for computing prediction coefficients as represented by a block 52. The 
prediction coefficients are based on past spectrum values. The prediction algorithm could be 

10 a standard algorithm for predicting a stationary time series from the finite past (such as the 
Trench algorithm discussed herein below). Hie filter coefficients computed by the prediction 
algorithm of block 52 are tested against three tap constraints represented by a block 50. If the 
constraints are not satisfied, an adjustment is made to R(0) by block 48 until the prediction 
coefficients do fall within the constraints. It may also be necessary to recompute the 

15 prediction coefficients with a constraint imposed. The prediction coefficients are inverted in 
polarity and used as the ICF tap weights for FIR filter 14 (FIG. IB). 

The R(0) adjust of block 48 can be used to establish an artificial noise floor for the 
algorithm. In such case, the R(0) adjust in essence simulates noise that either augments or 
reduces the actual noise on the communication channel and the algorithm responds 

20 accordingly to calculate ICF tap coefficients that have a built in positive or negative margin. 

An algorithm shown in FIG. 11 carries out the constraints represented by block 52. 
The autocorrelation function produced by autocorrelation block 46 is processed by a step 
represented by a block 60 using the predictive algorithm of block 50 to compute an initial set 
of ICF tap weights to be tested by the constraints. The ICF tap weights resulting from block 

25 60 are compared with a constraint A at query block 62. In one embodiment, constraint A is 
that the lowest order tap magnitude be smaller than 0.25. 

If constraint A is satisfied, i.e., passes, then the operation proceeds to a query block 
64, where the process determines whether constraints B and C are satisfied. In one 
embodiment, constraint B is that the sum of the tap magnitudes is smaller than 3, and 

30 constraint C is that the absolute values of the real and imaginary parts of all prediction taps (L 
and Q) is less than 2. 

If both constraints B and C are also satisfied, then the ICF taps are output and the 
algorithm is finished. If one or both are not satisfied, i.e., fail, however, then operation 

15 
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proceeds to a step 65, and the R(0) value of the autocorrelation function is adjusted. This 
adjustment changes the input autocorrelation function. Operation then proceeds back to the 
step represented by block 60 to recompute the prediction coefficients and repeat the described 
steps. 

5 However, if at query block 62 constraint A is not satisfied, then the operation 

proceeds to a step 66 and a constraint flag is then set, which from then on results in the 
invocation of a constrained algorithm, which incorporates the magnitude of the first tap as a 
constraint on the solution. Operation then proceeds to a query block 68, where the process 
determines whether constraints B and C are satisfied. If so, then the operation proceeds back 
10 to the step represented by block 60 to recalculate the prediction coefficients using the 
constraint. 

If at query block 68 one or both of constraints B and C are not satisfied, then 
operation proceeds to a step 65, the R(0) value is adjusted, and operation then proceeds back 
to step 60. 

15 In one embodiment of the invention, an iteration limit is enforced on the number of 

times the prediction coefficients are computed pursuant to block 60. For example, the 
prediction coefficients may be computed up to five times in an effort to satisfy the three 
constraints. If after the selected number of iterations, the constraints are still not satisfied, 
then the previous ICF tap weights are allowed to persist. 

20 In a farther improvement, rather than beginning the design iteration cycle of FIG 1 1 

afresh with each updated spectrum estimate, a means of incorporating information from 
previous designs is employed to shorten time required to calculate the prediction coefficients. 
This is shown in FIG. 12, in which the results of the initial tap weight computation, together 
with the results of the initial computation and last computation of the previous design, are 

25 used to determine the parameters for the second computation of the current design. The 
previous design refers to the final ICF taps output from block 64 after the last calculation and 
the current design refers to the ICF taps being derived in the course of the present calculation. 

In FIG. 12, the elements used in the algorithm of FIG. 11 are assigned the same 
reference numerals. In one embodiment, a state memory 70 stores the following state 

30 variables from the previous design (filter setting generation operation): 

• Whether constraint A failed on 1 st iteration of previous design. 

• Whether constraint A failed on any iteration of previous design. 

• First R(0) adjustment of previous design. 
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• Cumulative R(0) adjustment predicted by previous design for use in current 
design. 

These state variables, or other selected suitable variables used in the algorithm, are 
taken into account in computing the R(0) adjustment of the current design, as represented by 
5 a double ended arrow 72 between blocks 70 and 65. In addition, a leakage factor is included 
in the memory characteristics of state memory 70, as represented by an arrow 74. The 
leakage factor determines the weighting given to the past values of the variable relative to the 
present values of the variables. If the leakage factor is zero, the past has no influence, i.e., it 
is ignored. If the leakage factor is unity, the influence from the past is fully weighted. 
1 0 Normally a value slightly less than unity is used for the leakage factor. As represented by a 
broken line arrow 76 and arrow 72 to state memory 70, the state variables are transmitted to 
state memory 70. 

In a further improvement of the algorithm of FIG. 1 1, as represented in FIG. 12, by a 
query block 78, the iterations continue as time permits, further refining the noise floor / R(0) 

1 5 value. The R(0) value may be further increased or reduced. For example, there may be time 
for three iterations of the algorithm of FIG. 1 1 before the arrival of a new spectrum estimate. 
The algorithm may initially adjust the noise floor on the first iteration by a large amount for 
maximum convergence speed. The second iteration will then meet the constraints with 
excess margin. This will lead to a lowering of the noise floor on the third iteration, which as 

20 a result may or may not satisfy the constraints. The next design can then benefit from the 
previous design. In this fashion, a design satisfying a new interference spectrum can be 
found and applied quickly. The most recent design that passed the constraints is the one to 
use in the ICF. Iterations that are more recent, but did not pass the comtraints, or did not 
complete due to lack of time, are not used. However, the most recent iteration, regardless of 

25 whether it passed the constraints, is stored in the state memory for use by the next design. 

First, an initial unconstrained computation is performed in the manner described in 
FIG. 11. The results of this initial tap weight computation, together with the results of the 
initial computation and last computation of the previous design, are used to determine the 
parameters for the second computation of the current design. For example, in one 

30 embodiment there may be time for three iterations of Figure 1 1 before the arrival of a new 
spectrum estimate. The algorithm may initially adjust the noise floor on the first iteration by 
a large amount for maximum convergence speed. The second iteration will then meet the 
constraints with excess margin. This will lead to a lowering of the noise floor on the third 
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iteration, which as a result may or may not satisfy the constraints. The next design can then 
benefit from the previous design. In this fashion, a design satisfying a new interference 
spectrum can be found and applied quickly. 

The unconstrained algorithm employed to compute the filter settings is described in a 
5 paper: William F. Trench, "Weighting Coefficients for the Prediction of Stationary Time 
Series From the Finite Past," SIAM Journal of Applied Math, Vol. 15, No. 6, November 
1967, pp 1502-1510 (Hereinafter 'Trench", 'Trench paper" and/or 'Trench operations"). 

In using the Trench operations, it is necessary to begin with the autocorrelation 
function of the process that is to be predicted (i.e., predicted and then canceled in our 
10 application) when designing prediction taps using the Trench operation equations. We 
develop an estimate of the desired autocorrelation function in our ICF processing by inverse 
Fourier transforming the power spectrum of the ingress (tracked, masked, and otherwise 
processed to be made most suitable for noise canceling). When N taps are to be used in the 
prediction (or canceling) filter, plus the initial or feed-through tap, then the autocorrelation 
1 5 function is needed for displacements of zero to N, i.e., R(n), n= 0 to N. In the traditional case 
of designing prediction taps using the equations from Trench, the recursion formulae 
presented in Trench's Theorem 1 suffice, and the N+l values of the autocorrelation function 
are all that is necessary as input to the recursions. 

In the case of the constrained first tap magnitude operations, however, the traditional 
20 design procedure is not sufficient. If one design of the prediction taps is completed in the 
traditional way, and it is found that the first tap magnitude is larger than desired (constraint 
A), for whatever reason, then a remedy for this situation is as follows. Assume that the first 
tap in the prediction filter just derived is given as pi, and that it is desired to limit its 
magnitude to pljimit. We then desire to find the best values for the remaining prediction 
25 taps, p2 to pn, given that the first tap now must be constrained. 

We call upon Theorem 2 in Trench to solve this more difficult problem than the 
original, unconstrained design. As in the traditional unconstrained design, we use our 
. autocorrelation function, R(n), n= 0 to N, as the input to the recursions of the reference. Our 
autocorrelation function R(n) is used to define the autocorrelation matrix ® r . s , for r,s = 0 to N- 
30 1, of Reference 1, in the traditional unconstrained design. 

To accommodate the constrained design, we make the following modifications: . 
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1 . In Equation 9 of the Trench paper, we now define the matrix O r . s , for r,s - 0 to 
N-2, using R(n), n= 0 to N-2. Further, in Equation 9 of Reference 1, we define r| r , for r= 0 to 
N-2, as 

ilk = R(k+2) - [R(k+1) pj/lp^KplJimit], k= 0 to N-2. 
5 2. Then, the recursions of the Trench paper, Theorem 2, are carried out. 

3. After completing the recursions of the Trench paper, the solution to Equation 
9 is developed, and is denoted in the reference as §om, £im, • • • Um, where we have m = N-2 in 
our constrained case here. 

4. The prediction tap coefficients for this constrained tap design case are then 
10 pi,con S tr a med=lpi/|pi|][pUimit]; and 

Pn,constrained = 4n-2,N-2j for tt= 2 to N. 

5. The N values p^constmined, n = 1 to N, constitute the N coefficients for the 
desired prediction filter. Note that since noise canceling is our desire, the additive inverse of 
these coefficients is used in the ingress canceling filter. These operations conclude one 

1 5 embodiment of the constrained algorithm. 

FIG. 8 is a logic diagram illustrating operations performed corresponding to the filter 
weight computation algorithm of FIG. 7B in a different manner. The operations of FIG. 8 are 
performed for each sampling interval. Operation commences by sampling the 
communication channel to produce the sampled signal (step 802). If the communication 

20 channel is sampled while a signal of interest (desired signal) is present, the signal of interest 
is removed from the sampled signal (step 804). The sampled signal is then converted to the 
frequency domain to produce a spectral characterization of the sampled signal (step 806). An 
example of such a spectral characterization is illustrated in FIG. 5. 

The spectral characterization of the sampled signal is then compared to a threshold to 

25 determine whether the sampled signal is valid (step 808). One particular technique for 
performing this determination was described with reference to block 30 of FIG. 7B. If the 
spectral characterization is discarded, the operation of FIG. 8 for the particular sampled 
signal is complete and operation returns to step 802 wherein the commumcation channel is 
sampled during the next sampling interval. 

30 If the spectral characterization is valid, one or more tracking functions are applied to 

each spectral component of the spectral characterization (step 810). Then, the user mask is 
applied to each spectral component (step 812). These operations were described in detail 
with reference to FIG. 7B. Next, the noise floor of the spectral characterization is raised to a 
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budgeted SNR (step 814). This operation was described with particular reference to FIGs. 5, 
6, and 7B. All of these operations of step 810, 812, and 814 are performed in the frequency 
domain. The result of these operations is the production of a modified spectral 
characterization. The modified spectral characterization is then converted to the time domain 
5 (step 816). 

With the modified spectral characterization in the frequency domain, a first set of 
filter coefficients for the sample interval is computed (step 818). This operation may be 
performed according to the Trench paper. With the filter coefficients determined, a 
determination is made as to whether constraint A is passed (step 820). In the described 

10 embodiment, constraint A is passed if the first tap coefficient magnitude (after monic tap) is 
less than 0.25. If constraint A is passed, the filter coefficients are compared to constraints B 
and C. In the described embodiment, constraint B considers whether the sum of the tap 
magnitudes is less than 3.0 and constraint C considers whether the max tap I or Q component 
is less than 2.0. If all of these constraints are met at steps 820 and 822, the filter taps are 

15 applied to FIR filter 14 of FIG. 7A. 

If, however, constraint A is not met, all subsequent filter coefficient determination is 
performed using the constrained algorithm described above herein (block 824) and operation 
proceeds to decision block 826. If constraint A is met but constraints B and/or C are not met, 
operation proceeds to decision block 826 but the unconstrained algorithm continues to be 

20 employed. If a time out condition has occurred at decision block 826, e.g., four iterations, 
actual time constraint, etc., operation ends without determination of new filter coefficients. 
In such case, a previous set of filter coefficients are employed. 

If the time out condition of decision block 826 has not occurred, the modified spectral 
characterization and/or filter tap settings produced may be modified using prior 

25 operations/prior solutions (block 828). Then, the noise floor is adjusted (block 830) and 
either the unmodified or the modified algorithm is employed to compute again the filter 
coefficients (block 818). The operations of FIG. 8 will continue to be performed for the 
current sample until either a solution is produced or a time out condition is met. 

FIG. 13 is a system diagram illustrating the manner in which the present invention 

30 may be employed in a cable telecommunication system servicing high-speed data 
communications. In the cable telecommunication system, a head end unit 1302 couples to a 
plurality of cable modems 1304A-1304E via respective cable media 1306A-1306E. The head 
end unit 1302 services high speed data communications for the plurality of cable modems 
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1304A-1306E by acting as a gateway to a data network, such as the Internet 1308, a 
telephone network 1312, or other services network 1314. The head end unit 1302 also acts to 
distribute cable TV signals via the cable media 1306A-1306E to TV set top boxes (not 
shown). These cable TV signals therefore share the cable media 1306A-1306E with the high- 
5 speed data communications also serviced by the head end unit 1302. 

The head end unit 1302 services the high-speed data communications within a 
frequency band of interest. Interfering signals may reside within this frequency band of 
interest. Once source of these interfering signals may be radio frequency (RF) coupled from 
an RF source 1308. HAM radio operators, for example, produce radio frequency emissions 

10 in the frequency band of interest that may couple to the cable media 1306A-1306E as 
interfering signals. RF coupled interfering signals are typically intermittent However, some 
interfering RF signals may be predicted. Further, the cable system itself may produce 
interfering signals via adjacent channel interference or via infrastructure components mat 
operate in a faulty manner. These signals may be intermittent or may be somewhat continual, 

1 5 depending upon the cause of the interfering signals. 

Thus, according to the present invention, the head end unit 1302 and/or the cable 
modems 1304A-1304E employ the filtering operations of the present invention to remove the 
interfering signals from the communication channel serviced by the cable media 1306A- 
1306E. By performing these operations, the communication channel will service data 

20 communications at a higher effective data rate, thus increasing the throughput supported by 
the cable telecommunication system. 

In general, the described algorithms could be implemented in software that operates 
on a special purpose or general-purpose computer or in hardware. If the calculations 
necessary to compute the filter tap coefficients have not been completed by the next time 

25 interval T, the time interval is simply skipped and the algorithms operate on the narrow band 
interference energy in the following time interval T. 
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CLAIMS: 

1 . A method for canceling narrow band interference on a communication channel 
CHARACTERIZED by: 

sampling the communication channel to produce a sampled signal; 
5 spectrally characterizing the sampled signal across a frequency band of interest to 

produce a spectral characterization of the sampled signal, wherein the spectral 
characterization does not include a signal of interest; 

modifying the spectral characterization of the sampled signal to produce a modified 
spectral characterization; 
1 0 generating filter settings based upon the modified spectral characterization; and 

filtering the communication channel using the filter settings when the signal of 
interest is present on the communication channel. 

2. The method of claim 1, wherein in modifying the spectral characterization of 
15 the sampled signal to produce the modified spectral characterization, the method is further 

CHARACTERIZED by raising a noise floor of the spectral characterization to meet a 
budgeted signal-to-noise level. 

3. The method of claim 1, wherein in modifying the spectral characterization of 
20 the sampled signal to produce the modified spectral characterization, the method is further 

CHARACTERIZED by modifying spectral components corresponding to an expected 
interfering signal. 

4. The method of claim 1, wherein in modifying the spectral characterization of 
25 the sampled signal to produce the modified spectral characterization, the method is further 

CHARACTERIZED by modifying spectral components corresponding to a plurality of 
expected interfering signals. 

5. The method of claim 1, wherein in modifying the spectral characterization of 
30 the sampled signal to produce the modified spectral characterization, the method is further 

CHARACTERIZED by using spectral components of prior sampled signals. 
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6. The method of claim 5, wherein in considering the spectral components of 
prior sampled signals in modifying the spectral characterization of the sampled signal to 
produce the modified spectral characterization, the method is farther CHARACTERIZED by 
weighting more heavily the presence of an interfering signal in a prior sampled signal than 

5 the absence of the interfering signal in the current sampled signal. 

7. The method of claim 1, wherein the method is further CHARACTERIZED by 
sampling the communication channel to produce the sampled signal when the signal of 
interest is not present on the communication channel. 

10 

8. The method of claim 1, wherein Ihe method is further CHARACTERIZED by: 
sampling the communication channel to produce the sampled signal while the signal 

of interest is present on the communication channel; and 

removing the signal of interest from the sampled signal. 

15 

9. The method of claim 1 , wherein the method is further CHARACTERIZED by 
modifying the spectral characterization of the sampled signal to produce a modified spectral 
characterization in the frequency domain. 

20 10. The method of claim 1, wherein the method is further CHARACTERIZED by 

modifying the spectral characterization of the sampled signal to produce a modified spectral 
characterization in the time domain. 

1 1 The method of claim 1 , wherein the method is further CHARACTERIZED by 
25 discarding the sampled signal when total spectral density of the spectral characterization of 
the sampled signal exceeds a threshold. 
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12 A digital filter for canceling narrow band interference on a communication 
channel CHARACTERIZED by: 

a filter coefficient computation block that operates to: 

sample the communication channel to produce a sampled signal; 

spectrally characterize the sampled signal across a frequency band of interest 
to produce a spectral characterization of the sampled signal, wherein the spectral 
characterization does not include a signal of interest; 

modifying the spectral characterization of the sampled signal to produce a 
modified spectral characterization; and 

generating filter settings based upon the modified spectral characterization; 

and 

a digital filter that filters the communication channel using the filter settings when the 
signal of interest is present on the communication channel. 
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